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Résumé. - L’empoissonnement des barrages par des mugilidés est le principal moyen adopté en Tunisie pour 
assurer la continuité de la pisciculture continentale. La détermination des périodes d’abondance de certaines 
espèces revêt une priorité pour la réussite des opérations de récolte des alevins de muges. En effet, la res¬ 
semblance entre les alevins de Mugilidae et la dépigmentation de certaines espèces posent un problème pour 
l’identification des spécimens. Les analyses phylogénétiques ont montré que les populations de muge dépig¬ 
mentées sont les plus proches génétiquement de l’espèce Liza aurata. Elles ont aussi montré la robustesse 
d’un critère de différenciation entre Liza ramada et L. aurata, fondé sur la présence d’une ou deux bandes 
transversales de chromatophores à la base des rayons de la nageoire caudale. L’analyse spatio-temporelle des 
alevins de muges pêchés dans les différentes stations montre que l’abondance d’une espèce donnée, favorable 
à la récolte des alevins, est limitée dans le temps et dans l’espace. Pour Mugil cephalus, la période favorable 
à la récolte s’étale de fin septembre à mi-décembre ; pour L. aurata elle est entre janvier et mars, alors qu’elle 
est comprise entre février et avril pour L. ramada. 


Abstract. - Identification, spatiotemporal distribution and genetic characterization of grey mullet fry in Tuni- 
sian dams. 
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Stocking dams by Mugilidae is the principal mean adopted to ensure continuity and development of fish 
farming in Tunisia. Knowledge of the periods of abundance of each species is essential for the successful 
capture of mullet fry and for avoiding confusion between species. The great similarity between the mugid fry 
and the depigmentation of some juvéniles are the main issues. Genetic analysis was practiced to obtain a mor- 
phological criterion to differentiate Liza ramada from L. aurata and to identify depigmented fry. The spatio¬ 
temporal analysis of mullet fry caught in the various stations showed that the station Raoued is characterized 
by a high abundance of Mugil cephalus in September, while L. aurata is dominant in October. From January 
to May, L. ramada is abundant with the presence of L. saliens only in January. At Medjerda, L. aurata was 
abundant in October. M. cephalus abundance was low and that of L. ramada high from February to April. For 
stations P4 and drawbridge, L. ramada is commonplace between March and April. In the stations Hesyene and 
Maftouh El Maki, L. ramada is abundant in April. The phylogenetic analyses showed that the depigmented 
juvéniles belong to L. aurata. They also showed the presence of two or single band of chromatophores at the 
beginning of the caudal fin is a robust criterion to discriminate between L. aurata and L. ramada. The analy¬ 
sis of mullet fry in different stations shows that fry abundance was limited to short periods. M. cephalus is 
abundant from late September to mid-December and from January to March for L. aurata , while it is between 
February and April for L. ramada. 


Les muges sont des poissons de la famille des mugilidés 
qui vivent en bancs dans les eaux côtières. La structuration 
de ce taxon, malgré les avancées permises par les travaux 
moléculaires, est encore mal résolue (Durand et al., 2012). 
Les nombreuses études réalisées sur les muges ont révélé 
une large capacité d’adaptation à des milieux très différents. 
Ce sont des poissons euryhalins et eurythermes qui fréquen¬ 
tent les lagunes et les embouchures d’oueds, où ils passent 
la majeure partie de leur vie trophique. La large réparti¬ 
tion géographique des muges les soumet à une activité de 


pêche intense utilisant des engins de capture variés. De plus, 
l’aquaculture ne cesse de se développer et constitue une acti¬ 
vité à intérêt socio-économique soutenu (Gautier et Husse- 
not, 2005). L’élevage des muges sur de telles bases avait déjà 
été proposé en Tunisie par Pillai (1975). Les Mugilidae sont 
représentés par 25 genres (Durand et al., 2012). Cinq espè¬ 
ces seulement existent en Tunisie : Mugil cephalus, Chelon 
labrosus, Liza aurata, L. saliens et L. ramada. 

Les mugilidés présentent un grand intérêt économique 
et, même si les quantités pêchées sont loin d’être complè- 
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tement répertoriées, les captures en Méditerranée pourraient 
atteindre plusieurs milliers de tonnes (FAO, 2010). En Tuni¬ 
sie, des efforts sont actuellement déployés pour développer 
les activités d’élevage en milieu continental, en particulier 
par l’empoissonnement des retenues de barrages par diffé¬ 
rentes espèces. Cette action constitue une base possible pour 
l’établissement d’une activité piscicole extensive, ce qui per¬ 
mettrait de créer des emplois dans les régions défavorisées et 
d’approvisionner les populations riveraines en poissons frais. 
Les muges sont parmi les principales espèces exploitées dans 
les barrages tunisiens qui constituent des écosystèmes assez 
proches du milieu naturel des alevins. Leur élevage repose 
sur la récolte des alevins en milieu naturel. Les espèces les 
plus recherchées pour l’aquaculture sont M. cephalus et 
L. ramada en raison de leur euryhalinité et de leur croissance 
rapide en eau douce. En revanche, L. aurata et C. labrosus 
ont une croissance moindre dans les retenues de barrages 
(Heldt, 1948). L. saliens est une espèce déconseillée en rai¬ 
son d’une croissance trop lente et d’une taille maximale peu 
importante (Sarig, 1981 ; Ardizzone et al., 1988 ; Cataudella 
et al., 1988). La récolte des alevins de muges se fait par cap¬ 
tage ou pêche des alevins, au niveau des communications 
mer-lagune, par des équipes de pêcheurs spécialisés du Cen¬ 
tre Technique d’Aquaculture (CTA) qui assure aussi l’ense¬ 
mencement de 9 millions d’alevins par an dans 25 grands 
barrages en plus des lacs collinaires. Les alevins capturés 
doivent être immédiatement identifiés, acclimatés et trans¬ 
portés directement vers les retenues de barrages pour être 
ensemencés. L’identification des alevins de muges doit être 
précise afin d’alimenter les barrages présentant des caren¬ 
ces en certaines espèces. Pour stimuler le développement de 
cette activité aquacole et donner plus d’envergure à l’aqua¬ 
culture continentale et pour que cette opération puisse durer 
(FAO, 1983), il faut résoudre certains problèmes d’identifi¬ 
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cation des espèces d’alevins de muges ensemencés dans les 
barrages tunisiens. 

Les critères d’identification anatomique et morphologi¬ 
que des alevins de muges sont semblables à ceux des adul¬ 
tes, mais d’utilisation plus délicate. De nombreux critères 
distinctifs chez les adultes (membrane adipeuse de l’œil, lon¬ 
gueur relative des nageoires pectorales, épaisseur et structu¬ 
re de la lèvre supérieure,...) ne sont pas encore bien établis. 
Les difficultés les plus communes sur le terrain sont l’inca¬ 
pacité d’identifier des spécimens totalement dépigmentés 
et l’absence de caractères morphologiques distinctifs entre 
L. aurata et L. ramada. De plus, l’identification de certaines 
espèces, notamment celles du genre Liza, basée sur la pig¬ 
mentation (Cambrony, 1984), n’a en fait jamais été vérifiée. 

L’objectif de ce travail est de fournir une clé d’identifica¬ 
tion des alevins de muges facilement utilisable sur le terrain. 
Pour cela, nous avons utilisé une analyse génétique pour dif¬ 
férencier les alevins de L. ramada de ceux de L. aurata, ainsi 
que pour l’identification des alevins dépigmentés issus de la 
station Pont Mobile (PM). De plus, la répartition spatio-tem¬ 
porelle des alevins de muges dans la région nord-est de la 
Tunisie a été analysée. 

MATÉRIEL ET MÉTHODES 

Les alevins de muges ont été collectés au niveau des 
embouchures de la majorité des oueds, des lagunes et des 
effluents côtiers dans la région nord de la Tunisie (Fig. 1). 
Le site de collecte le plus important est situé au niveau de 
l’oued Medjerda, il regroupe la station “Pont Mobile” 
(PM) (37°01’08 N, 10°04’31 E) et celle de l’embouchure 
de l’oued (37°00’42 N, 10 o ll’34 E). Les principaux autres 
sites de collecte sont : Maftouh El Makki (37°03’20 N, 
10°10’23 E), porte 4 (P4) (37°04’25 N, 10°09’52 E) et les 


Figure 1. - Carte des zones de collectes 
des alevins de muges. [Map of collec¬ 
tion areas for the mulletfry.] 
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Tableau I. - Données relatives aux séquences utilisées lors des analyses phylogénétiques. A : Océan Atlantique ; M : Mer Méditerranée ; P : 
Océan Pacifique. [Sequence data used inphylogenetic analyses. A: Atlantic Océan; M: Mediterranean Sea; P: Pacific Océan.] 


Espèces 

Code Genbank 

Site de collecte 

Références 

Liza aurata 

EU122432 

EU122433 

EU122434 

M : Tunisie,îles Kuriat 

Trabelsi et al., 2008 

EU224057 

A: France, Bretagne Sud 

www.fishtrace .org 

EF439541 

M: Espagne 

www.fishtrace .org 

Chelon labrosus 

EF427545 

EF427544 

M: Espagne, mer Cantabrique 

www.fishtrace .org 

DQ197935 

A : Espagne, îles Canaries 

Pescabase: Genetic catalog of Spanish marine fish 

Liza saliens 

EU715507 

EU715505 

A : Espagne 

Heras et al., 2009 

Liza ramada 

EU224058 

EU224059 

A : France, Pertuis d’Antioche 

www.fishtrace .org 

Mugil cephalus 

EU122430 

M : France, golfe du Lyon, Sète 

Trabelsi et al., 2008 

EU036449 

EU036450 

M : Grèce, Kavala 

www.fishtrace .org 

EU036450 

P : Chili 

Fraga et al., 2007 


stations de Raoued (36°59’05 N, 10°12’22 E) et oued Hesie- 
ne (37°02’28 N, 10°09’23 E). Les alevins ont été capturés à 
l’aide de filets du type senne italienne modifiée ou “trattina”. 
Deux techniciens à pied maintiennent les extrémités du filet, 
encerclant le banc d’alevins préalablement repéré a l’œil nu. 
Ce dernier est dirigé vers le sac du filet où il sera accumulé 
pour assurer un transfert rapide. Les alevins sont alors immé¬ 
diatement récoltés et déversés dans une cuve de 1 m 3 , bien 
oxygénée. La quantité d’oxygène injectée est dépendante de 
la charge en alevins. L’eau utilisée pour le transport provient 
toujours de la station de collecte des alevins. Ces derniers 
sont transportés jusqu’aux barrages où ils seront déversés. 

Une fois arrivés au barrage, les alevins sont acclimatés 
par le déversement de quelques seaux d’eau en provenance 
du barrage. Cette méthode simple a montré son efficacité 
pour l’acclimatation des alevins de muges. Après 30 min 
permettant de minimiser les chocs thermique et osmotique 
en équilibrant la salinité et la température, les alevins sont 
déversés dans le barrage. Un échantillon représentatif de la 
quantité totale d’alevins capturés (généralement 1%) est col¬ 
lecté dans des bouteilles qui sont congelées (-20°) en atten¬ 
dant leur identification morphologique et moléculaire au 
laboratoire. 

Analyse moléculaire et génétique 

L’ADN est extrait en utilisant la méthode de Trabelsi et 
al. (2002). Sur chaque individu conservé dans l’alcool, nous 
avons prélevé environ 0,15 cm 2 de la nageoire caudale. Une 
région de 376 paires de base du gène codant le cytochrome 
b, a été amplifiée par PCR en utilisant les amorces : Sens 
5’-AGCCTACGAAAACCCACCC-3 ’ (Meyer et al., 1990) 


et Antisens 5 ’ -AAACTGCAGCCCCTCAGAATGATATT- 
TGTCCTCA-3’ (Cantatore étal., 1994). 

Pour chaque individu, les réactions de PCR sont effec¬ 
tuées dans 50 pA de mélange réactionnel contenant 50 ng 
d’ADN matrice, 0,2 pM de chaque amorce, 2,0 U. HiTaq 
polymérase, 0,2 mM de dNTPs et 5 pl de tampon (Bioprobe). 
L’ADN amplifié est vérifié par migration sur un gel d’aga- 
rose, puis séquencé. Les séquences sont alignées à l’aide de 
BioEdit (Hall, 1999). Les analyses de reconstruction phylo¬ 
génétique sont effectuées par la méthode du maximum de 
vraisemblance (Saitou et Nei, 1987) avec une matrice basée 
sur la distance à deux paramètres (Kimura, 1980) présente 
dans le logiciel de construction d’arbres phylogénétiques 
MEGA 4.0.2 et moyennant une approche par recherche heu¬ 
ristique de PAUP (Swofford, 1998). 

Les analyses de la reconstruction phylogénétique sont 
effectuées en se basant sur plusieurs séquences du gène du 
cytochrome h extraites de la banque de gènes appartenant 
aux espèces suivantes : M. cephalus, L. ramada, L. aurata, 
L. saliens et C. labrosus (Tab. I). 

Des séquences de cytochrome b issues des populations 
de muges des stations PM et P4 ont été obtenues pour les 
groupes d’alevins suivants (isolats) : 

- Isolat 1 et 2 PM-DPTUN 2011 (alevins dépigmentés de 
la station PM) ; 

- Isolat 1B-PM TUN 2011 et Isolat 2 1B-P4 TUN 
2011 (alevins présentant une seule bande de chromatophores 
au départ de la nageoire caudale) ; 

- Isolat 1 2B-PM TUN 2011 et Isolat 2 2B-P4 TUN 
2011 (alevins présentant deux bandes de chromatophores au 
départ de la nageoire caudale). 


Cybium 2013,37(1-2) 


69 
















Identification et répartition spatio-temporelle des alevins de muges 


Mili et al. 


Analyse spatio-temporelle 

L’Analyse en composantes principales (ACP) a été appli¬ 
quée uniquement aux stations Raoued et Medjerda qui pré¬ 
sentent une fréquence de présence d’alevins de muges lar¬ 
gement répartie dans le temps. De plus, ces deux stations 
présentent le stock le plus important en alevins de mugilidés 
dans la région nord de la Tunisie. Ainsi, les échantillons col¬ 
lectés ont été suffisants pour réaliser des analyses statistiques 
fiables. Le reste des stations de collecte d’alevins n’ont pas 
fait l’objet d’échantillonnages suffisants pour permettre de 
déterminer avec précision les fréquences de présence, en rai¬ 
son de leur présence ponctuelle et très limitée dans le temps. 

Pour déterminer de manière précise les périodes d’abon¬ 
dance pour chaque espèce, on a procédé à une Analyse en 
composantes principales (ACP) en utili¬ 
sant le logiciel STATISTICA 8.0. 


RESULTATS 

Analyse génétique 

Le complexe des muges est divisé 
en trois clades ((BP) > 99%). Si l’on ne 
tient pas compte des alevins, deux clades 
ne comprennent que les séquences d’une 
seule espèce (M. cephalus et L. ramada ), 
et le dernier clade se subdivise en trois 
groupes monophylétiques constitués par 
les espèces L. saliens, L. aurata et C. 
labrosus (Fig. 2). Chaque groupe monos¬ 
pécifique est très fortement soutenu statis¬ 
tiquement (BP = 100). Le deuxième clade 
de l’arbre phylogénétique est uniquement 
constitué de séquences de L. ramada de 
France et des haplotypes appartenant aux 
populations présentant une seule bande 
au départ de la nageoire caudale (stations 
Porte 4 et Pont mobile). Ces populations 
forment un groupe monophylétique qui 
présente des pourcentages de bootstrap 
élevés. Le dernier groupe est divisé en 
deux sous-groupes dont l’un contient 
les séquences de populations de muges 
dépigmentées provenant de la station PM 
et L. aurata d’Espagne qui présente des 
valeurs de bootstrap de l’ordre de 96%. 
Le deuxième sous-groupe est constitué 
d’haplotypes appartenant aux muges pré¬ 
sentant une barre transversale de chro- 
matophores, doublée par une barre en arc 
de cercle, des stations P4 et Pont mobile, 
ainsi que L. aurata de Tunisie et de Fran¬ 
ce. La présence d’une barre transversale 


de chromatophores doublée par une barre en arc de cercle au 
départ de la nageoire caudale pour L. aurata, et la présence 
d’une seule bande transversale de chromatophores au départ 
des rayons de la nageoire pour L. ramada constituent donc 
un critère de différenciation des deux espèces. Les haplo¬ 
types des alevins se groupent avec les séquences des deux 
espèces du genre Liza ne modifiant pas la valeur respective 
de bootstrap (BP : 100). Il n’y a aucune différenciation signi¬ 
ficative des populations de différents sites de collecte ce qui 
peut être dû à des flux de gènes entre les populations. Pour 
les populations dépigmentées et celle de L. aurata , l’absence 
de bifurcation selon les analyses du maximum de vraisem¬ 
blance permet clairement de regrouper ces deux populations 
(BP = 96%). 
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Figure 2. - Arbre phylogénétique en utilisant la méthode du maximum de vraisemblan¬ 
ce. (Les valeurs de bootstrap (BP) sont données pour chaque nœud si la valeur est supé¬ 
rieure à 50 en pourcentage). [Phylogenetic tree using the maximum likelihood method. 
Bootstrap values (BP) are providedfor each node if the value is greater than 50 per¬ 
cent).] 
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Pour l’analyse phylogéographique des alevins de muges, 
les séquences partielles de cytochrome b de populations tuni¬ 
siennes présentant une seule bande au départ de la nageoire 
caudale se positionnent avec celles de l’espèce L. ramada. 
Les populations présentant une barre transversale de chro- 
matophores doublée par une barre en arc de cercle ainsi que 
les populations dépigmentées se positionnent avec l’espèce 
L. aurata. Une petite portion d’un gène mitochondrial appro¬ 
prié peut donc permettre de positionner spécifiquement des 
alevins de muges. 

Analyse spatio-temporelle des alevins de muges 

Les deux composantes principales contribuent pour 
79,83% de la variance totale de l’ACP pour les poissons 
de la station Raoued (Figs 3,4). La première composante 
principale, qui contribue pour 51,41% est caractérisée par 
une contribution positive et élevée de L. ramada et par une 
contribution négative et élevée de M. cephalus. La deuxième 
composante principale, qui contribue pour 28,42% dans la 
variance totale, est caractérisée par une contribution négative 
et faible de L. aurata et de L. saliens. Le mois de septembre 
est caractérisé par une abondance de l’espèce M. cephalus. 
En octobre, L. aurata est l’espèce abondante (Figs 3,4). Les 
mois de février, mars, avril, mai, novembre et janvier sont 
caractérisés par une abondance de l’espèce L. ramada alors 
qu’en janvier l’espèce L. saliens présente une abondance 
relativement importante. 

A la station Medjerda, les deux composantes principales 
de l’ACP contribuent pour 100% de la variance totale (Figs 
5,6). La première composante principale, qui contribue pour 
99,99% de la variance totale, est caractérisée par une contri¬ 
bution positive élevée de L. ramada et négative et élevée 
de L. aurata. Selon cette analyse, les mois de février, mars 
et avril peuvent être groupés ensemble. Le mois d’octobre 
est caractérisé par une grande abondance de L. aurata et par 
la présence de M. cephalus en faible quantité. Pendant la 
période allant de février à avril, L. ramada est l’espèce la 
plus présente. 

DISCUSSION 

Nos analyses moléculaires ont permis de mettre en place 
une clé d’identification facile à appliquer sur le terrain par 
les responsables des opérations de pêche d’alevins avant 
de passer à l’ensemencement dans les retenues de barrages. 
Nos résultats sont cohérents avec plusieurs autres études qui, 
sur la base de paramètres morphologiques et moléculaires, 
ont permis d’identifier des juvéniles de muges (Trabelsi et 
al., 2008). La structuration phylogénétique est comparable 
à celle qui est présentée dans les autres publications sur la 
phylogénie des mugilidés (Aurelle et al., 2008 ; Heras et al., 
2009 ; Durand et al., 2012). Toutefois, de nombreuses ques- 
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Figure 3. - Contribution des espèces étudiées dans la Composante 
principale appliquée à la station Raoued. [Loading ofspecies on the 
Principal Components applied to Raoued station.] 



Figure 4. - Contribution des mois dans la Composante principale 
appliquée à la station Raoued. [Loading ofmonths on the Principal 
Components applied to Raoued station.] 


fions concernant le positionnement intra-générique et intra- 
spécifique restent sans réponse et seules des études portant 
sur un nombre élevé d’individus d’origine différente et com¬ 
binant l’utilisation de marqueurs nucléaires (microsatelli¬ 
tes) peuvent permettre de résoudre les positionnements aux 
niveaux spécifiques voire populationnels de ce taxon. Des 
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Figure 5. - Contribution des espèces étudiées dans la Composante 
principale appliquée à la station Medjerda. [Loading ofspecies on 
the Principal Components applied to Medjerda station.] 



Figure 6. - Contribution des mois dans la Composante principale 
appliquée à la station Medjerda. [Loading ofmonths on the Princi¬ 
pal Components applied to Medjerda station.] 


travaux de ce type ont par exemple permis l’identification de 
trois espèces biologiques au sein d’individus du complexe 
M. cephalus à Taiwan (Shen et al., 2011). Les marqueurs 
moléculaires peuvent permettre de détecter des isolements 
reproductifs au sein d’espèces cryptiques vivant en sympa- 
trie. Cela permet aussi d’infirmer ou de confirmer des intro- 
gressions suggérées par l’analyse mitochondriale. 
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Les conséquences des observations sur l’exploitation des 
alevins dans la station de Medjerda et Raoued sont de plu¬ 
sieurs ordres. Tout d’abord, il apparaît que les abondances 
favorables à la récolte des alevins d’une espèce donnée sont 
souvent limitées à des périodes restreintes. Pour M. cepha¬ 
lus, la période favorable à la récolte s’étend de fin septembre 
à mi-décembre, ce qui correspond à la fin de la période de 
reproduction de cette espèce (Farrugio, 1975 ; Brusle, 1981 ; 
Brusle et Brusle, 1977). Dans cette zone d’étude, la période 
d’abondance maximale des alevins de muges coïncide avec 
celle qui a été observée précédemment par Vidy et Franc 
(1987) et par Chauvet (1986). La diminution de la période 
d’abondance au cours de l’année étudiée (2011) peut s’ex¬ 
pliquer par une surexploitation des ressources en alevins ou 
encore par les changements climatiques ou la pollution dans 
cette zone (Ennouri, 2012). Pour L. aurata, le pic de collecte 
se situe entre janvier et mars, alors qu’il s’étale entre novem¬ 
bre et avril dans les études précédentes (Chauvet, 1986 ; 
Vidy et Franc 1987). La période entre février et avril corres¬ 
pond au pic de collecte de L. ramada dans les eaux tunisien¬ 
nes alors que les études antérieures le signale entre janvier et 
avril (Farrugio, 1975 ; Chauvet, 1986 ; Vidy et Franc, 1987). 
L. saliens n’a pas présenté d’abondances exceptionnel¬ 
les dans la station de Medjerda mais elle présente de fortes 
abondances en septembre pour la station de Raoued. Les tra¬ 
vaux signalant des abondances importantes de cette espèce 
indiquent que la meilleure période de pêche s’étale du mois 
de mai au mois d’octobre. Pour la station P4, L. ramada est 
l’espèce la plus abondante dans les échantillons prélevés 
entre mars et avril. En ce qui concerne les stations Maf- 
touh El Maki et Hesyene, L. ramada est abondante en avril, 
alors qu’elle est présente au niveau du pont mobile entre le 
mois de mars et avril.La détermination précise des périodes 
d’abondance de chaque espèce est primordiale pour la réus¬ 
site des opérations de récolte et pour éviter la confusion des 
espèces. 

En conclusion, l’identification moléculaire des divers 
membres de la famille des Mugilidae présente un grand inté¬ 
rêt fondamental mais aussi appliqué, car ces poissons sont 
parmi les principales espèces exploitées dans les barrages 
tunisiens. Ces espèces représentent la majeure partie des ale¬ 
vins ensemencés dans ces retenues d’eau et leurs débarque¬ 
ments sont les plus importants. La présente étude a permis 
de mettre en place une clé d’identification, facile à appli¬ 
quer sur le terrain en prenant en considération la répartition 
spatio-temporelle et les critères morphologiques simples 
des alevins de chaque espèce. Ce travail a fourni une base 
de données pour l’identification morphologique des alevins 
de muges en Tunisie. La documentation des périodes et des 
lieux, où les alevins de chaque espèce sont les plus abon¬ 
dants, revêt une importance primordiale pour la réussite des 
opérations de récolte. 
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